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Li et al. reply: In their Comment on our recent 
Letter 0, Bhattacharya et al. note that for the geom- 
etry used in their experiment || to try to detect the 
modes supported by the Meissner state of a d-wave su- 
perconductor predicted by Zutic and Vails £9, a field of 
first vortex penetration H* cl of 300 Oe was measured. 
They then argue that the crossover field H* separating 
the regime dominated by nonlocal electrodynamics and a 
local regime is of order only 20 Oe, suggesting that non- 
local effects can be neglected for their geometry for most 
of the intermediate field range H* < H < H^. Instead, 
they attribute their failure to observe these modes to the 
absence of gap nodes due to an out-of-phase admixture of 
a secondary order parameter component. Here we ques- 
tion whether the theory of Zutic and Vails can be applied 
in this intermediate range, argue that a d + is state in 
YBa 2 Cu3 06.95 (YBCO) is unsupported by other types 
of experimental evidence, and continue to maintain that 
the nonlinear Meissner effect is unlikely to be observed in 
any of its manifestations, due in part to nonlocal effects. 

In our paper, we pointed out simply that the nonlocal 
crossover field H* ~ ^ (Ao is the penetration depth 

for supercurrents in the plane, and £o and £oc are the co- 
herent lengths in the plane and along the c axis, respec- 
tively), for the geometry in question, with a dominant 
(001) surface and H || a, b, is generically of the same order 
of magnitude as the field below which the Meissner state 
is thermodynamically stable, H c \ ~ 4 7r ^ °x where 
Aoc is the penetration depth for supercurrents along the 
c axis, and R — ^kokq c with kq and kq c the plane and the 
c-axis Ginzburg-Landau parameters, respectively. Bhat- 
tacharya et al. argue (l]|| that the experimentally mea- 
sured field of first vortex penetration is as much as an 
order of magnitude higher. With this problem in mind, 
we in fact alluded to precisely this difficulty of determin- 
ing the correct lower critical field in our paper Q . 

We do not have a complete understanding of the high 
value of H* ± consistent with all experimental observations 
at this time. An obvious explanation is the existence of a 
surface barrier, as is commonly observed in type-II ma- 
terials ||. On the other hand, the expected hysteresis 
in such systems due to barrier profile asymmetry has ap- 
parently not been observed. 

The real questions are as follows: For fields H& < 
H < is the sample still in the pure Meissner state 
and can one expect to observe manifestations of the non- 
linear Meissner effect? We believe not. Although the 
bulk magnetization has not occurred, the intermediate 
field state is thermodynamically unstable and vortices 
may be easily trapped at the sample corners or within 
the skin depth. Two recent high-resolution penetration 
depth measurements |||7j failed to measure the predicted 
temperature dependence in this field range, but did ob- 
serve a large linear-i/ field dependence most likely at- 
tributed to trapped vortices || . These vortices may make 
the field penetration layer highly nonuniform, leading to 
a broadening of the Zutic- Vails harmonics. 


Bhattacharya et al. argue jjJD that their failure to 
observe the predicted small nonlinearities is strong evi- 
dence for the existence of a true gap in the bulk YBCO 
system, as created, e.g. by a bulk d + is order parameter 
with Ag/Ad of order a few percent. This claim is however 
inconsistent with several other experiments which indi- 
cate the existence of a residual normal fluid in YBCO, 
m one case down to 30mK @. Were a small out-of-phase 
order parameter component A s , of order several degrees 
K, to exist, the gap in the density of states would ex- 
ponentially supress quasiparticle excitations at mK tem- 
peratures. On the other hand, a much smaller gap, of 
order tens or even hundreds of milliK is insufficient to 
eliminate the signal in the experiment of Ref. || . 

The experiment of Bhattacharya et al. || is a clever 
and careful attempt to measure a very small effect, and, 
as pointed out by Zutic and Vails |4|, in principle, is a 
sensitive probe of gap structure. However, the interpre- 
tation of the results is complicated by the influence of 
trapped flux, and it seems prudent that effects of this 
type be ruled out or understood before claiming the ob- 
servation of more exotic effects, such as the d + is state 
invoked by Bhattacharya et al. 
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